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We proposed the spontaneous CP-violation in the Simplest Little Higgs model. In 
this model, the pseudoscalar field can acquire a nonzero vacuum expected value. It 
leads to a mixing between the two scalars with different CP-charge, which means 
spontaneous CP-violation happens. It is also a connection between composite Higgs 
mechanism and CP-violation. Facing the experimental constraints, the model is 
still alive for both scenarios in which the extra scalar appears below or around the 
electro-weak scale. We also discussed the future collider tests on CP-violation in 
the scalar sector through measuring h2ZZ and hih2Z’ vertices (see the definitions 
of the particles in the text) which provides new motivations on future ete~ and pp 
colliders. It also shows the importance of the vector-vector-scalar- and vector-scalar- 


scalar-type vertices to discover CP-violation effects in the scalar sector. 
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I. INTRODUCTION 


The discovery of a 125 GeV Higgs boson [I| 2] by the ATLAS and CMS collaborations 
[3] in 2012 implies the success of the standard model (SM) because the measured signal 
strengths are consistent with those predicted by the SM [4] B|. However, the electro-weak 
symmetry breaking (EWSB) mechanism is an important topic and researches on physics 
beyond the SM (BSM) are still necessary and attractive. 

For example, to solve the little hierarchy problem, Arkani-Hamed et al. proposed the 
Little Higgs (LH) framework [6] in which the collective symmetry breaking (CSB) mechanism 
[6] was used to forbid the quadratic divergences in the Higgs potential at one-loop level. 
The LH framework contains a lot of models. All of them are special kinds of composite 
Higgs models [7] thus each of them must contain a global symmetry which is spontaneously 
broken at a high scale f > v where v = 246 GeV is the vacuum expected value (VEV) of 
the Higgs field. The SM-like Higgs boson is treated as a pseudo-Nambu-Goldstone boson 
corresponding to one of the broken generators and EWSB is generated dominantly through 
quantum correction thus the Higgs boson can be naturally light [6} [7]. Usually the gauge 
group is also enlarged thus there are extra gauge bosons with their masses at O(f) scale. 
LH models are effective field theories (EFT) below a cut off scale A ~ 4r f. Below the 
scale A, a LH model is weakly coupled, but we do not know what would happen above 
A. Among those models, the simplest Little Higgs (SLH) model [SHIO] has the minimal 
extended scalar sector in which their are only two scalars. In the SLH model, a global 
symmetry [SU(3) x U(1)]° is spontaneously broken to [SU(2) x U(1)]? at scale f. The gauge 
symmetry is enlarged to SU(3) x U(1) and spontaneously broken to the electro-weak (EW) 
gauge symmetry SU(2); x U(1) at scale f as well. And at the EW scale v, the gauge 
symmetry is further broken to U(1)em like what happens in the SM. If CP-violation is 
absent in the scalar sector, one of the scalars is the SM-like Higgs boson (denoted as h), and 
the other is a pseudoscalar. 

CP-violation is another important topic in both SM and BSM physics. In 1964, CP- 
violation was first discovered through the Kz — mr rare decay process [11]. More CP- 
violation effects have been discovered in K- and B-meson sectors [2]. All these measured CP- 
violation effects can be successfully explained by the Kobayashi-Maskawa (K-M) mechanism 


[12] which was proposed by Kobayashi and Maskawa in 1973. They showed that a nontrivial 


CP-phase can appear in the quark mixing matrix (named CKM matrix [13]) if there 
exist three generations of fermions. However, the succeed of K-M mechanism is not the end 
of CP-violation studies. For example, the observed matter-antimatter asymmetry in the 
universe [2] [14] requires new sources of CP-violation because SM itself cannot generate such 
a large asymmetry [6]. Thus it is attractive to study new CP-violation sources. Till 
now, the scalar sector is still an unfamiliar world for us and there may be lots of hidden 
new physics, including new sources of CP-violation. Thus in this paper, we focus on extra 
CP-violation in the scalar sector. 

Theoretically, there are already many extensions of the SM which contains new CP- 
violation sources. For example, if we add more complex scalar singlets or doublets, there 
may be CP-violation in scalar sector which can leads to a CP-mixing Higgs boson 
1. Some of these models may be CP-conserving at the Lagrangian level and CP-violation 
can arise only from a complex vacuum, which was called the spontaneous CP-violation 
mechanism [19]. This mechanism was proposed by Lee in 1973 [19] as the first kind of two- 
Higgs-doublet model (2HDM) [17]. Moreover, spontaneous CP-violation mechanism is also 
a possible solution of the strong-CP problem [23], and it may have further connection with 
lightness of the Higgs boson as well [24]. Besides these models, spontaneous CP-violation 
in the scalar sector can also arise from the composite framework. There are already two 
examples, one is the next-to-minimal composite Higgs model (SO(6)/SO(5), or equivalently 
SU(4)/Sp(4)) [25], and the other is the Littlest Higgs model (SU(5)/SO(5)) [26]. In each 
model, CP-violation occurs when the pseudoscalar field acquires a nonzero VEV. In this 
paper, we will propose the possibility of spontaneous CP-violation in the SLH model through 
the realization of the same mechanism. This model can also appear as one of the candidates 
to solve strong-CP problem as mentioned above. More details on this topic will appear in a 
forthcoming paper [27]. 

Phenomenologically, we can test new CP-violation effects directly or indirectly. The 
indirect effects may appear in the electric dipole moments (EDM) of electron and neutron 


[28], modifications in meson mixing parameters [29], or anomalous ZZZ couplings [80]; while 


1 For the 125 GeV Higgs boson, LHC measurements preferred a CP-even one and excluded a CP-odd one 
at over 30 level through the final distribution of h > ZZ* — 4@ decay assuming no CP-violation in the 
Higgs interactions |22|. However, a CP-mixing Higgs boson is still allowed since the contribution from 
pseudoscalar component should be loop suppressed. 


the direct effects may be discovered in ht*r~ or htt vertices through measuring the final 
state distributions [31]. If another scalar is discovered and we denote the scalars as hi» (hi 
is the SM-like Higgs boson and hə is the extra scalar), we can also discover CP-violation 
in the scalar sector through directly measuring tree-level vector-vector-scalar- (VV S-) and 
vector-scalar-scalar- (V S S-) type vertices, such as hyaVV and Vhyhg vertices 7, according 
to the CP-properties analysis [24]. Based on this idea, the author and his collaborators 
recently proposed a model-independent method to measure the CP-violation effects in the 
scalar sector through etet — Z* — Zhy,Zh2,h hz associated production processes at 
future ete~ colliders [32]. In that research, the product of the three vertices was used as a 
quantity to measure the magnitude of CP-violation [33]. However, in the SLH model, 
the author and his collaborators recently showed the Zh h2 vertex is suppressed by a factor 
(v/f which means it is difficult to test. Thus to test CP-violation in the SLH model, 
we can turn to extra heavy gauge bosons for help. 

As a summary, the model studied in this paper is attractive both theoretically and phe- 
nomenologically. This paper is organized as following: in we briefly review the 
CP-conserving SLH model, build the SLH model with spontaneous CP-violation, and obtain 
the domain interactions; in [section TIT] we consider the constraints on this model, especially 
in the scalar sector; in [section IV] we discuss the tests on CP-violation effects in this model 
at future ete” or pp colliders; and in section V] we present our conclusions and further dis- 
cussions. In the appendix we also presented the improved SLH formalism [34] 
which is very helpful for the model building. 


II. MODEL CONSTRUCTION 


In this section, we first briefly review the CP-conserving SLH model and then construct 
the spontaneous CP-violation SLH model. We will also derive the useful vertices in the 
spontaneous CP-violation SLH model. In both models, we have the same nonlinear realiza- 
tion for Goldstone bosons. We also have the same particle spectra in both models, while 
in the CP-violation model, the scalars are both CP-mixing states. The CSB mechanism 


and loop corrections in the Higgs potential are also similar in both models. The difference 


2 Here V denotes a massive gauge boson. For the SM gauge group, V = W or Z; while for LH gauge 
groups, V can also denotes extra heavy gauge bosons. 


comes from an extra explicit global [SU(3) x U(1)]? breaking term which is absent in the 


CP-conserving model. 


A. A Brief Review of the CP-conserving SLH Model 


The SLH model contains two scalar triplets ®;5 which transform as (3,1) and (1,3) 
respectively under the global [SU(3) x U(1)]? transformation [SHIO] [35]. At a scale f >v, 
[SU(3) x U(1)?? breaks to [SU(2) x U(1)}? and ten Nambu-Goldstone bosons are generated, 
eight of which should be eaten by massive gauge bosons during spontaneous gauge symmetry 
breaking SU(3) x U(1) > SU(2), x U(1) > U(1)em. Two physical scalars are finally left. 


The nonlinear realized scalar triples can be written as 


ð = ei pits 01x2 ; ® = e e7i9/ts 01x2 . (1) 
feg fsg 
where p is a mixing-angle between the two scalar triplets. The matrix fields © and ©’ are 
separately 
Q= 1 į bxs ae , and O= 1 (oes + ney (2) 
f\ v2 g 0 F\ v2 gt 0 


in which ¢ = ((v, +h —iG)/v2, ga is the usual Higgs doublet and y = (y?, £7)" is 
another complex doublet for Goldstones corresponding to heavy gauge bosons following the 
conventions in [85]. 


The covariant derivative term is 


where 
Dab- (4) 
g is the weak coupling constant and the gauge fields matrix is [8] [9] 
1 wt y® 
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=|} -1 |+ 1 +— | wW- X- 
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where Oy is the EW mixing angle 3 and complex fields Y°(Y°) = (Y'+iY?)/V2. The 
terms including G,,G“ in (3) give the masses of gauge bosons. Before EWSB, v, = 0; while 
after EWSB, vp is generated through quantum correction. It must be close to v and their 


difference arises at O ((v/f)”) level. To the leading order of (v/f), we have [35] 
gf 


my = > and mx = my = ==. (6) 


V2 


The other three neutral degrees of freedom will mix with each other at leading order of (v/f) 


through the matrix [35] 


A —Sw SxCw CxCw A’ 
Z = Cw SxSw CxSw A® (7) 
Ze 0 Cx —Sx B 


ji 


where 0 = arcsin (tw/v3). The corresponding masses at leading order of (v/f) are then 


gu 2 
=0 = d 1=4] z 7g. 8 
MA » mz Dey? and mz 3 Tp (8) 


The massless gauge boson is photon. If we go beyond leading order of (v/f), the gauge 


bosons will have further mixing with each other. For example, in charged sector, WF and 


X+ will mix with each other at O((v/f)°) level, and W(X)* will acquire their relative 


mass corrections at O ((v/f)”) level. While in neutral sector, the off-diagonal elements of 


the mass matrix Mj, in the basis (Z, Z’, Y?) are nonzero. Using an orthogonal matrix R, it 


can be diagonalized as (RMj,R*) 


— ? 
= MpÔpq Where m, are the gauge bosons’ masses. The 


pq 
neutral gauge bosons acquire their mass corrections as 


bm? mZ, 4 m2 6 

d N, @ , se (>) out 0 @ (9) 
mz f Mz f My2 f 

We denote the corresponding mass eigenstates as Z, Z’, and Y. Their mixing 


2 3 
Rzz ~O (6) ) 5 and Rzy2 N Raye ~O (G) ) . (10) 


A and Yt do not participate further mixing. A mass splitting between Y' and Y? is also 


generated as ôm? /m? ~ O ((w/ fÊ). 


3 In this paper, we denote Sa = sina, Ca = cosa, and ta = tana, for any angle a. 


The six neutral scalar degrees of freedom can be divided into CP-even (h and yt) and CP- 
odd (n, G, G’, and y?) parts where y°(g°) = (yt + iy?) /V2. A straightforward calculation 


showed that after EWSB, the kinetic terms can be written as 


1 
Lrin = z (Ə!hð,h + "yt y" + K,j0"G;0,,G;) (11) 


with G; runs over the four CP-odd scalar degrees of freedom and K;; 4 ôi; means the CP- 
odd part is not canonically-normalized *. To find out the canonically-normalized basis, we 
should consider the gauge fixing term together. The two-point transitions between gauge 


bosons and scalars arise from the cross terms of 0,0; and G,®;. These transitions can 


be parameterized as VE Fpi G; and their contributions should be canceled by CA er ee 
from the gauge fixing term. It can be checked straightforwardly that (see appendix|section A] 


or for more details) a new basis 


(ii, Gp) a : n a) (12) 


is canonically-normalized. Gp is just the corresponding Goldstone of Vp. 


In the fermion sector, each left-handed doublet must be extended to a triplet thus there 
must be additional heavy fermions. In lepton sector, a heavy neutrino N; should be added 
for each generation. While in the quark sector, choosing the “anomaly-free embedding” [B6], 
T with Q = 2/3 is added as the parter of t, D and S with Q = —1/3 are added as the parters 
of d and s separately. The Yukawa interactions are then 


Ly = i Np lL, — A ag det (1, 2, Lx) 
y (arahat it Atak »®}) i= i dry det (©, 5, Q3) 
H (indin? + Aada?) Qn Eang det (8,050); (13) 
where the left-handed triplets are 


L= (vr, br, iN)” Qı = (dz, —uz,iDz)’ , 


i ? 


Qo = (51, —cL i91), Q = (tr, br, iT)" . (14) 


The first line is for leptons where lr; runs over (e, u, T)r; the second line is for the third 


generation of quarks where dpr; runs over (d, s,b, D, S)r, and the last line is for the first 


4 Details on the improved formalism to treat this case can be found in the appendix [section A]and [34]. 


two generations of quarks where ur, runs over (u,c,t,T)r. A ~ 4rf is a cut-off scale. A 
right-handed quark with index a or b must be a mixing state between an additional quark 
and its SM partner, for example, i are mixing states between tr and Tgr. To the leading 


order of (v/f), The heavy fermions’ masses are [9] 


mi =X fse, mo = y [Asca + [dbsal?F, (15) 


for Q = T, D, S and q =t,d(d,), s(d2). To the leading order, the corresponding partners in 

SM sector have the masses 

of wt |AZA2| qv 
"V2 flages + [assaf Y2 


CSB mechanism keeps all neutrinos massless °. Other fermions require their masses (simi- 


(16) 


larly, to the leading order) 


v U 


aa: 
= ; Mp = A ; Mu, = U,C> 
A Jon” ° 42r K ? Ia 


in which y% are eigenvalues of matrix A% and Yuc are eigenvalues of matrix AÎ¥. To this 


J 
My 


(17) 


step, we ignored small mixing between q and Q. Consider this kind of mixing Ajg, a mass 
correction dm,/m, ~ O (A?9/M8) is generated. 

Last, let’s turn to the scalar potential. In the discussions above, we assume the Higgs 
doublet acquire a correct VEV to derive the particle spectra everywhere. However, at tree- 
level, GEA term is forbidden due to the CSB mechanism. The Higgs potential can be 


generated through Coleman-Weinberg mechanism [87] at loop-level as 
SVa = —ôm? (hth) + 6 (hth). (18) 


The CSB mechanism forbids quadratic divergence in thus [8] 


3 K Pm, X g?m3, (ltt), A 
( m Tio 8n? ( pip 2 mZ 1 Ym 8 2 me, y (19) 
(5m), oop 3 4 mp l 
— l TaD. a) 
gt i me 1 F+) 1 my 1 2 
2 (ym 2) _ oe) (20) 


5 In the first term of (13), we can also use ®; instead of 62, but we cannot have both terms together if we 
assume massless neutrinos. If we perform this replacement, m4, in should also be changed to A% fcp. 


Here A ~ 4rf is a cut-off scale and A; = J2m: /v, which means the contributions from the 
first and second generations of fermions are ignorable. When mr is heavy enough, EWSB 
can be generated through these loop corrections. 
Now the pseudoscalar ņ is still massless due to an accidental global U(1) symmetry. 
Adding a term 
6V = =pl, + Hc. (21) 


in the potential °, 7 acquires its mass 


2 2 


m2 = — cos ( Z ) pak 22 
7 speg V2fseca/ Sate ve 
and the Higgs potential acquires another correction 
(Va), = — (8m?) (hth) + (6d), (hh) 
sm e 23 
m ( ) a Pee ( ) 3 ( ) 


Two-loop contributions to 6m? can be absorbed into the possible contributions from un- 
known physics at the cut-off scale A [6}[38] which can be parameterized as (5m). toop = —cf?. 
We can roughly estimate |c| ~ O(107?). 


B. Spontaneous CP-violation in the SLH Model 


In (21), u-term provides the 7 mass. In general, u? can be complex, but its argument 
can always be absorbed into the shift of 7 (which is equivalent to a rotation of ®;). Besides 
this, 7 cannot acquire a nonzero VEV, thus there is no CP-violation in the scalar potential. 
Comparing with the CP-conserving case in we can add another term and 
becomes 


2 
SV = old, + € (eie) +e (24) 


Here € is also required to be small (for example, € < O ((v/f)?) thus the CSB mechanism is 


not significantly broken). In general, u? and € can be complex, but we can shift 7 to make at 


6 This term breaks the CSB mechanism explicitly which means a quadratic divergence in the Higgs potential 
can be generated at one-loop level. Thus numerically u should be very small comparing with f. In the 
convention of this paper (which is the same as that in [85]), the degrees of freedom in ©’ cancels with 
each other thus 7 dose not acquire additional mixing with y?. 
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least one of them real. If we choose u? real, when c€ is still complex, CP-symmetry would be 
explicitly broken in the scalar sector. However, if both u? and e are real, 7 is also possible 
to acquire a nonzero VEV which means spontaneous CP-violation happens. In this paper, 


we focus on the spontaneous CP-violation case. 


According to (24), denote a = vp/(W2fsacs), we have 


n A 22 2 v2n 
V, = — u’ f?sgcgCa COS (=) + ef °83C3C, COS (25) 
" = V2f sace fsace 


Minimize this potential, we found that when 
1? < def? |secacal ; (26) 
(7) = 0 becomes unstable thus 7 would acquire a nonzero VEV 


2 
u 
Un = = +V 2f sgcg arccos | —_-——— ], 27 
by = (1) = VB seq arccos ( zE) (27 
which means spontaneous CP-violation is possible. For simplify, we choose “+” in the 
equation above from now on. We denote € = v,/ (/2f SgCg), and the scalar mass term is 


M? ef s2q82¢ h 


Lm = -5 (hun) (28) 


Ef?’ 82qS2¢ de f’ Cs? n 
Here M should be close to 125 GeV and it includes all the quantum correction effects from 


and T. Nonzero off-diagonal elements means the mass eigenstates cannot be CP 


eigenstates. Define the mass eigenstates (in which hı is SM-like) 


h ca —S h 
oe 0 o (29) 
hy So Co n 
we have the mixing angle 
1 2e f? S2452 
d= arctan (z = de ft (30) 
and scalar masses 
M? + 4e f? 5? M? — 4e f? 82 
M2 = y J £ + ( 5) E cog + cP san 8350 š (31) 


T These quantum corrections are not affected by the CP properties of the scalar sector which means 
and derived in the CP-conserving model can be simply transported into the CP-violation case 


11 


We can see that only when both u? and € are nonzero, CP-violation can occur, which means 
in this model, CP-symmetry is also collectively broken ° 

For the Yukawa couplings, we can also choose all the couplings real thus there is no explicit 
CP-violation. Complex CKM matrix can arise from the mixing between a SM quark and an 


extra quark, which is the same mechanism as that in {18}. 


C. Some Useful Interactions in this Model 


In the CP-violation SLH model, mixing between h and 7 can modify some of the vertices 


in the CP-conserving model. The hVV couplings can be parameterized as 


gv ms 
Livy = Ea > (Evh + č vhə) vv") (32) 


where V denote the mass eigenstates. For real vector fields, V* = V. To the leading order 


of (v/f), we have 


P 7 c 7 i s 
Č,w =C Cw=S6, Č,z = -Čz = T C27 = =C = aa (33) 
Cw 
286 U 
Č -], 34 
e ae AE J. ax = Spe (>) (34) 
2 Qcecow 6) 7 2secow e) 
C1, Y = — j C2 Y = A (35) 
3t3 Cy f Er! f 


Here č; x comes from the X* and W= mixing [85], while ¢;y comes from Z, Z’, and yY? 
mixing. They arise at O ((uv/f)°) which are extremely small comparing with &; w/z/z'. 


For the antisymmetric type Vhn couplings °, we parameterize it as 


Ly nana = $ (hdl he — hah) (Cran Zu + Cinna Zp + Franaa) (36) 


The results to the leading order of (v/f) are 


1 3 24/2 
Cohths = (z) , Co hiha = v2 (>) , CY hiha = a (37) 
2/2c3,top \ f V3 —teytos \Í 


8 The case € absents was already discussed above. The case u? absents allows a nonzero vy, but € = 7/2 
that the off-diagonal elements in are still zero. A shift of ņ (rotation of ©) can remove this € hence 
it is trivial. A nontrivial € requires nontrivial u? and e. 

9 We don’t consider the symmetric type couplings (h10“h2 + h20“h1) here because they cannot contribute 
anything in the processes with in-shell gauge boson(s). 
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which are the same as the CP-conserving case, since h10”ho — h2d"h, = hô”n — no”h. 


The scalar trilinear interactions should be 
i! 2 1 2 
Ls = — 5 At2at ihe = graf hahi; (38) 


where to the leading order of (v/f), the dimensionless coefficients 


2ES2a (3C2¢ — 1 2€S8¢ (3C2a — 1 
nn = a (1 — 353) alae) 5, (2 — 953) Ba 
528 S28 
Qec? À 
S26 f 
2 3 a 1 2 Q 3 — 1 
A211 = Cp (1 = 353) vi2es2¢ (3c2 ) — sọ (2 = 382) V 2€ 820 (3¢2¢ — 1) 
= 598 = 
x ‘ se 
+6c2s9—— + bys? V282 (40) 
f S28 
A in the equations is the Higgs self-coupling constant. 
m The Yukawa couplings for SM leptons and quarks f = £, q can be parameterized as 
a m Z 
= Ly = >. = ((c1,phi + C2, pho) fifr) + H.c. (41) 


f 


For f = u,c,b,v,£, the pseudoscalar degree of freedom dose not couple to these fermions, 


thus we have 


ar =C and p= So; (42) 


while for q = d, s,t, the coupling coefficients 


U C28 + C26p U C28 + C26p 


f V 2596 I V 2898 l 


Here 6, = —1 for the third generation (q = t) and 6, = +1 for the first two generations 


Ciq = Co | 10980 and C24 = $9 — 109Co (43) 


(q = d,s). The imaginary parts are generated by the left-handed mixing between light and 
heavy quarks. 0r = arctan (ae / d2) at the leading order of v/f is the right-handed mixing 
angle. Here we don’t consider the possible flavor changing couplings. The Yukawa couplings 
including a heavy quark should be 

fy =-->). a ( (c1.qhi + c2Qħ2) QLQR 
Q 


+q ((ciLqgħi + Cat qħ2) Pr + (cing + Cor qh2) Pr) Q + H.c.) (44) 
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where Pr/p = (14 7°)/2. The coefficients 

C19 = cor (=) Piga = sp (=) line. as 
Here 6g = +1 for the third generation (Q = T) and dg = —1 for the first two generations 
(Q = D, S), which different with those for SM fermions. s29, X dgm,q/mg thus for the first 
two generations, we have so, < 1. The other four coefficients including both light and 


heavy quarks are 


V (C28 — C20p)$20p , Sp Mq C28 — C29p . Maf 
c =g + iĝos = —ĵÎoCo — isg——, (46 
1L,q "OF p Q ! Jasa Q V2mo 528 mov ( ) 

V (C28 — C29pR)S29R . Sp Mq C28 — CR . Maf 
c = 5 — idoc = —098¢ + icg ; (47 
2L,q °F oe Q ore Q V2mo 598 mov ( ) 

2 
C28 + C9 . V C28 + C26p 
= ô — —1], 48 
ari ey aa; (( s25 ) ) 8) 
a2 a cap + Catm -1 (49) 
2R,q Q20 Vasa | oof Sa : 


In the calculation of cir, the improved formalism affects on their imaginary parts since the 
7 component in G cannot be ignored due to the improved SLH formalism [34]. For the third 
generation, m:/mr ~ O(v/f), thus cz, can reach O(1). But for the first two generations, 


Ma/Mo K v/f means Cipg < 1. 


II. RECENT CONSTRAINTS ON THE MODEL 


As a BSM model, SLH always face many direct and indirect constraints, such as collider 
searches for new particles predicted by the model and EW precision tests. The scalar sector 
contains an extra scalar hs, whose properties are quite different from the SM-like scalar. If 
it is light enough (mz < m,/2), it should also face the hı cascade decay constraint. As a 
model with new CP-violation source, we should also discuss the EDM constraints [28]. In 


this paper, we don’t discuss more details about quark flavor physics. 


A. Direct and Indirect Constraints on f 


In the SLH model, the modifications on S and T parameters are sensitive to the new 


scale f. Thus before LHC Run II, the S and T parameter constraint 41| on f used to 
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be the strictest one. f > (4—7) TeV at 95% C.L. when tg ~ (1 — 10) [42] [43]. In the SLH 
with spontaneous CP-violation, this constraint is similar, because the S and T parameters 
are note sensitive to mg and c2 w/z when cp w/z < 1. 


However, since LHC Run II began, the lower limits on exotic particles increase quickly 


hence the corresponding new physics scales are pushed higher. In the SLH model, X* and 
YY?) gauge bosons couple to SM fermions with a suppression factor v/f, thus they are 
difficult to be produced at LHC. However, couplings between Z’ and SM fermions have the 
same order with those in SM !°, thus Z’ searches at LHC can provide a direct constraint on 
f. Recently, using 36.1 fb~' luminosity at \/s = 13 TeV, ATLAS collaboration set a new 
constraint mz = 4.5 TeV at 95% C.L. [44] for the sequential standard model (SSM) in 
which Z’ couples to SM fermions with the strengths in the SM. 
In the SLH model, the signal strength 


(oz Brzse+e-) ory ~ 0.36 Kd/u + 1.14 
(oz Brz_se+e- SSM l Kd/u + 0.78 


u ~ 0.49, (50) 


in which 
—— J dx dao fa(x1) fq(@2)d(a122 — m32,,/s) 
d/u = 
/ f dxydxofu(v1) fa(t2)d(a1v2 — m2,,/s) 
for mz = (4—4.5) TeV, using the MSTW2008 PDF [46]. Comparing with the results shown 


in [44], it can be roughly estimated f = 7.5 TeV at 95% C.L. 


~ (0.2 — 0.25) (51) 


Comparing with the indirect constraints discussed above, we can see that the Z’ direct 
searching experiments can provide the strictest constraint on f in the SLH model for most 
p region. Only in the case with extremely large (such as 2 10) or small (such as < 0.1) tg, 


the indirect constraints can be stricter. 


B. Constraints on the Properties of Extra Scalar hə 


hy couples to SM particles dominantly through its h component, since the couplings 
between 7 component and SM sector are highly suppressed by the high scale f. Experimen- 
tally, for a light h2, it mainly face the direct searches through ete~ + Zhə at LEP; while 
for a heavy he, it mainly faces the direct searches through gg > ha > W*W7/ZZ at LHC. 


Both production cross sections a suppressed by a factor s3. When mə < mı/2, it should 


10 These couplings are the same in CP-conserving and CP-violation models. 
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also face the hı — 2h2 rare decay constraint. Theoretically, the allowed parameter region 
also depend on the details of EWSB. 
For mz ~ (15—80) GeV, experimentally, LEP direct searches through ete~ > Z* + Zh 


associated production process gave [47| 
s9 < (0.1 — 0.2) (52) 


at 95% C.L. assuming Bra, = 1. €zn,n. cannot be constrained at LEP since it is sup- 
pressed by a factor (v/f)?. When mz < m,/2, it must face the hı rare decay constraint as 
well. In SLH model, the dominant exotic decay channel is hı — 2h2 with a branching ratio 


Brah = Vnys2n,/T1. The partial decay width 


XT” 4mž 
Phish = F i= Pre (53) 


while the hı total decay width 
Py = Phr + V'1sM- (54) 


Based on the Higgs signal strengths measurements using full 2016 data [4] B|, we perform a 


global-fit and obtain an estimation 
Brexo L 0.2 and sọ < 0.4, (55) 


at 95% C.L., which is a bit stricter than the previous constraint from LHC Run I [48}. 
We show the branching ratio distribution in According to the figures, when 
mz ~ (20 — 60) GeV, we have sg < (0.04 — 0.16) which is a stricter constraint than that 
from LEP direct searches. The numerical results are not sensitive to f and £. 
Theoretically, the allowed parameter region also depends on the details of EWSB, es- 


pecially the contributions from cut-off scale, 6m? = —cf*. In the CP-violation case, (23) 


becomes 


nth)”, (56) 


e (2È — 2coac 
ÖV} = 2e f? (2cacrace — CtoCre) (nth) 4 ( : 7 i x) ( 
383C% 
leaving the other contributions to ôV, unchanged. For f = 8 TeV, in the light h2 scenario, c 
is favored in the region (0.01 — 0.02) since larger cz is excluded by the Higgs data. However, 


if c < 0.01, EWSB requires larger sọ which was excluded by the Higgs rare decay constraints, 
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FIG. 1: The hy —> 2h2 decay branching ratio distribution in sg-mz plane with f = 8 TeV. From 


left to right, we choose tg = 1,3,6, respectively. 
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thus smaller cp would lead to the exclusion of light ho scenario. Larger f requires smaller c, 
for example, if f = 12 TeV, the lower limit of c reaches about 4 x 1073. 

For a heavy hə (with ma 2 200 GeV), experimentally it is constrained by LHC direct 
searches. At LHC, the gluon fusion process acquires dominant contribution through top 
quark loop, and the amplitudes through heavy quark loops are suppressed by (v/f)?, so 


thus orn /On sM © $3. If mg < 2m, the branching ratios of hy are the same as those of a 


SM-like Higgs boson with the mass mg. For mz > 2m, another decay channel hy — 2h, 


d2 f? Am2 
Tr — 2U fy 5, 57 
hə—>2hı 32TMoə m3 ( ) 


Its branching ratio can reach (20 — 30)% when mz = 300 GeV. If mz Z 350 GeV, ho > tt 


opens with a partial width 


decay channel can also open. Recently, ATLAS collaboration performed the direct searches 
through the channels pp > ha > W+W-—, ZZ for mz > 200 GeV with 36.1 fb! luminosity 
at ,/s = 13 TeV [49] 50|. If m < 1 TeV, the strictest constraints come from the hy > ZZ 
decay channel. Comparing with the SM theoretical predictions [52], we have a rough 


estimation 
(0.1 — 0.4), for me ~ (0.2 — 0.3) TeV; 


so Š 4 0.2, for m2 ~ (0.3 — 0.7) TeV; (58) 
(0.2 — 0.4), for mə ~ (0.7 — 1) TeV. 
at 95% C.L. These constraints are a bit weaker than the light mz region. 
Theoretical constraints here are similar to those in the case with light hg. c ~ (0.005 — 


0.03) is favored in the heavy hz scenario. In this scenario, the results are not sensitive to 
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f or 8. Bound on mz is sensitive to c, but not sensitive to sg, which is different from the 


properties in light hz scenario. 


C. EDM Constraints 


The EDM effective interaction can be written as 
idẹ -+ 
LEDM = = (59) 


which violated P- and CP-symmetries. In the SM, CP-violation comes only from complex 
CKM matrix so that the leading contributions to the EDMs of electron and neutron arise 


at four- and three-loop level respectively. It is estimated that 
desm ~ 1078 e-cm, — dngm ~ 10° e- cm, (60) 
both of which are far below the recent experimental constraints [53] [54] 
Ide| < 8.7 x 107” e- cm, Ida] < 3.0 x 107” e- cm, (61) 


at 90% C.L. However, in some BSM models, electron or neutron EDM can be generated at 
one- or two- loop level, which means it may face strict experimental constraints. 

In the SLH model with spontaneous CP-violation, the leading contribution to electron 
EDM comes from the two-loop “Barr-Zee” type diagrams with F = t,7T,D,S running 
in the loop, see the left diagram in Following the calculations in [56], we have 


FIG. 2: Dominant Feynman diagrams contributing to EDM of electron and quarks, and CEDM of 


quarks. F running in the loop includes t, T, D, S. 


Y Y g 
e,q eq d d q q 
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the analytical expression for the EDM of an electron as 


de 3G FAemMeSgCo C28 F C20R F mip mip 
= —— i = pna 2 
e (27)? (5 D3 D Qr $9 28 g m? g m2 ’ (6 ) 


in which the function 


g(z) = J ana — — In a a (63) 


Numerical results showed that de is not sensitive to the masses of extra heavy quarks. For 


0.2 < tg X 8, in the whole mass region mz ~ (20 — 600) GeV, we have 


8 TeV 
Ide] <8 x 107” ( Ja a) e-cma ft. (64) 


The constraints from electron EDM are not strict due to the suppressions by @ and f. 
For a neutron, its EDM comes from not only quarks’ EDM, but also their color EDM 
(CEDM) operator [28] [55] [56] 


igs 7 _ v a a 
OcepM = Fay" V(t Gia (65) 


where dy is the CEDM of the quark, t° denotes the color SU(3) generator, and i, j are color 
indices. The u quark EDM comes only from the left diagram in just like that 
for electron; while the d quark EDM acquire contributions from both the left and middle 
diagrams in because of the left-handed mixing between d and D quarks. The 
CEDM of quarks come from the right diagram in Calculate at the EW scale, the 
quarks’ EDM and CEDM in the SLH model with spontaneous CP-violation are [56] 


2m 
du D — de; 
3m, (66) 


Ma are v må m? 
t= gnc ape (P) CGC) 
a 
nt OE mie)! (a) (6+ 
+(e) Gn) = (0-2 
my me, 


7 Gra; ae C28 + rr p 
da 2a 


a-m- OCN) “ 
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in which the function 


1 
z 1 —2z(1— zx), z(1-— z7x) 
=- d l . 70 
al T a-r- Zz o 
After the running to hadron scale, the neutron EDM 
dn d 4 da : 
E E aay = 0.16— + 0.19du. (71) 
e 


Numerically, for 0.2 S ts < 8, in the whole mass region m2 ~ (20 — 600) GeV, we have 


8 TeV 
dal S 1.4 x 107” ( ; =) e-cm, (72) 


which is still below the experimental limit. The constraint from neutron EDM is weaker 
than that from electron EDM. 


Besides the “Barr-Zee” type diagram, there are also one-loop diagrams and Weinberg 


operator contributing to neutron EDM, see the Feynman diagrams in Following 


FIG. 3: Additional Feynman diagrams contributing to neutron EDM. 


g 
Y g 
A rane Le 


(46)-(49), we can estimate the one-loop contribution to neutron EDM (the left and middle 


diagrams in|Figure 3) as 


(73) 


0.7s9c9 Ma |M? — m2 m? 
|Sdn1-1oop| ~ = a - 2 l ( Bp , 


327tg = ufm?, u? 
where the scale 4 ~ O(v). This result is sensitive to mp and mg. For mp ~ O(f) and 
mz S O(v), |ôdn, 1100p) £ O07? — 10-7”) e-em. The Weinberg operator (the right diagram 
in Figure 3) (7) 


W rabera v,b up, 
Ow = = è Ge GG" (74) 


in which f*” is the structure constant of SU(3) group, contribute to neutron EDM as 


odnW (9.8 MeV)w. (75) 
e 
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In the SLH model with spontaneous CP-violation, we have 


7 Gras v CR t + C28 m? m? 
um Gar ee (we) a) = 


where the function [56] 


eo (1 — v)(we)? 
WN | i | eae lv) aaa 7) 


Typical |dd,w| < O(10-*8) e- cm. Thus we can conclude that for neutron EDM, the 
contributions from are sub-dominant. 


There are also upper limits on heavy atoms’ EDM. The recent measurement on !°Hg 
atom’s EDM set new limit |dy,| < 7.4 x 107% e- cm at 95% C.L. which provides an 
indirect constraint dn < 1.6 x 107% e- cm [59]. The SLH model with spontaneous CP- 


violation is still allowed by this new indirect constraints. The theoretical estimation on the 
EDM of Hg contains rather large uncertainties [60], thus it cannot directly provide further 


constraint on this model. 


IV. FUTURE COLLIDER TESTS OF THE CP-VIOLATION EFFECTS 


Recent Higgs data have already confirmed the 0+ component of hı [22]. Following the 
idea in [24} B2], we should try to measure tree level hyVV and hihgV vertices to confirm 


CP-violation in the scalar sector. For different hy mass, we need different future colliders. 


A. Measuring həVV Vertex 


If hy is light (for example, mo < v), it is difficult to be discovered at LHC because of 
the large QCD backgrounds at low mass region. To test this scenario, we need future ete 
colliders. For example, at CEPC or TLEP with \/s ~ (240 — 250) GeV, the hoa VV 
vertex can be measured through the ete~ — Z* — Zhə associated production process. Its 
cross section is [47| 


as Taz, (8st — 48%, +1) - så F m m? = 12m3 z me m? (78) 
? 96s(1 — m? /s)}2sfy ciy s’ s s s’ s ? 


in which the function 


F(x,y) = V1 +a? +y? — 2z — 2y — Qay. (79) 
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With 5 ab™! luminosity at CEPC, the inclusive discovery potential on sg can reach 5o if 
sọ ~ 0.15 at low mass region (mə < 70 GeV) through the “recoil mass” technique 
(64) [65]. This result does not depend on the decay channel of hz, and it is not sensitive 
to mz in this region. With the the help of “pr balance cut” method to reduce large 
backgrounds with photons, the 50 discovery bound on sg can reach about 0.1 with a tiny 
breaking of inclusiveness. If we completely give up the inclusiveness in this measurement and 
consider only the hy — bb decay channel, the 50 discovery bound on sg can be suppressed 
to about (4—5) x 107? according to 11. This result means the allowed regions obtained 
in are still possible to be discovered at 5a level at CEPC with 5 ab‘ luminosity. 
For larger mə when it is close to Z-peak, large ZZ background will decrease the sensitivity 
on sg measured though this channel. 

For light hz, we can also measure sg through Z > Z*(f f)h2 rare decay, if an e*e~ collider 
runs at Z-pole (ys = mz). The branching ratio [2] 


2 T mz=m2 3 202 a2 
s2 , q’ po _ mab sin“ 6 
Brz-+z*ho = H) sinoao | “We — m2 F (2 j {= B2 ) ; (80) 
Z 


where q is the invariant mass of Z*. The momentum of hg in initial Z frame and the relative 


velocity between h and Z* are respectively 


n C ae ) a 


MZP2 
j= , 82 
p3 + \/(p3 + m3) (p3 + @) 


With 10!? Z boson events as the goal of a “Tera-Z” factory, the typical sensitivity to this 
rare decay branching ratio is about (1078—1077) [67], which means it has a better sensitivity 
to discover nonzero sg comparing with the Zh associated production channel in the whole 
mass region m2 S< 70 GeV. 

For a heavy hə (for example, mz ~ O(v)), LHC future direct searches will discover it or 
set a stricter limit on sg, through its ZZ decay channel [68]. Through merely visible leptonic 
decay channel, with 3 ab~', the 5a discovery bounds would be around sg ~ (0.1— 0.2), which 
is similar to the current upper limits using the combination of hg > 4@ and hy —> 22v 


channels [49] [52] [68]. We also expect the 242v channel can help to increase the sensitivity 


11 Simulation details about the cross sections of the background channels were not shown in the text of [32]. 
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on sg at future LHC. When mə = 0.6 TeV, the 242v channel would become more sensitive 


than the 4£ channel [49]. 


B. Measuring hjh2V Vertex 


Based on the improved formalism of SLH model [84], we obtain the Zhh vertex in (37). 
CH, h 18 Suppressed by a factor (v/f)? < O(10~°) that the associated production channels 
cannot be used to measure this vertex. Similarly, precision measurements on hy > Z® hs 
are also useless to test this vertex, since the typical 50 discovery bounds for such rare decay 
channels are of O(1073) [69]. That means we must turn to the heavy neutral gauge 
boson sector for help. 

According to (37), CZinihy 1S Suppressed by a factor (v/f), and there is no suppression in 
CYhihj- These vertices will become helpful to confirm the 0~ component in at least one of 
the scalars. Since mz >> m1, the decay branching ratio 


m, v \? 
Brz = 2 83 
TAS Beg af ( itn) ee) 


Assuming the heavy quark masses mp > mz /2 thus Z' = FF decay channels cannot be 
opened. The total width Tz; + 6.5 x 10~3f if we choose the “anomaly free” embedding [9]. 


Numerically, we have 


8 TeV\? 
Diving airo E eer 84 
12 
ftg 


When 8 ~ 7/4, this decay channel vanishes, while if 8 is close to 0 or 7/2, there is an 
enhancement by ton It decreases quickly when f increases. 

For this process, we need future pp colliders with larger \/s, for example, (50 — 100) TeV 
[70]. Since at LHC, when mz = 5 TeV, the event number of pp > Z! + hıh cannot 
reach O(1) with 3 ab~! luminosity [70]. However, with the same luminosity at \/s = 100 TeV 
pp collider, the events number can reach Npp+z/+hjho ~ O(10? — 10%) for mz: ~ 5 TeV, and 
Nop+Z!-+haho ~ O(10 — 10?) for mz ~ 10 TeV [70]. This implies the Z'hıho vertex in the 
SLH model is testable at \/s = 100 TeV pp collider. 


If we can discover both nonzero vertices həaZZ and Z’hih2, we can confirm the CP- 


violation effects in the scalar sector. 
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V. CONCLUSIONS AND DISCUSSIONS 


We proposed the possibility of spontaneous CP-violation in the scalar sector of the SLH 
model in this paper. Through adding a new interaction term, € (oje) +H.c., in the scalar 
potential, the pseudoscalar field 7 can acquire a nonzero VEV which means CP-violation 
happens spontaneously. Both scalars then become CP-mixing states. In this paper, we 
denote hı as the SM-like Higgs boson with its mass mı = 125 GeV, and hə is the extra 
scalar. Based on the improved SLH formalism (see in the appendix), we derived 
the interactions in this model. 

Facing strict experimental constraints, the spontaneous CP-violation SLH model is still 
not excluded. LHC Run II data have already push the lower limit of the scale f to about 
7.5 TeV, which means the EW precision tests only provide sub-dominant constraints on f. 
For the extra scalar h2, we have two scenarios based on its mass, mz ~ O(v) or ma < v. 
For a light ha, the most strict constraint comes from hı — 2h, rare decay channel. The 
95% C.L. upper limit on sg is (0.04 — 0.16) for mə ~ (20 — 60) GeV. While for large 
Mz ~ O(v), the 95% C.L. upper limit on sg varies in the region (0.1 — 0.4), especially when 
mə ~ (300 — 700) GeV, the 95% C.L. upper limit on sọ is about 0.2. In both scenarios, tiny 
but nonzero contributions from cut-off scale are necessary. As a CP-violation model, it must 
also face the EDM constraints. Since the effects are suppressed by sgu/f, the constraints 
are weak. The most strict EDM constraint comes from electron, which favors 0.2 S tg S 8 
in the whole ma ~ (20 — 600) GeV mass region. 

We also discussed the future collider tests of this model. The basic idea is to discover 
nonzero hyaVV and Vhyhz vertices. For a light hə , we can test hyZZ vertex at future 
ete” colliders, as Higgs factories or Z-factory. With 5 ab~' at CEPC, for mz < 70 GeV, 
sọ ~ (4—5) x 107? can be discovered at 50 level; while with 10!? Z-boson events at 
Z-pole, we can have a better sensitivity in the same mass region. For a heavy hg with 
M ~ O(v), the vertex can be tested through gg > hy > ZZ channel at LHC. With 
3 ab™" luminosity, the 5a discovery bound is around (0.1 — 0.2) through merely the 40 decay 
channel. The 242v decay channel is also expected to help increase the sensitivity on sọ, 
especially in large mz region. Based on the improved formalism, we know the Zh,h2 vertex 
is suppressed by (v/ f)’, thus we must ask a heavy gauge boson, such as Z’, for help. Since 


Brzhh £ O(10~4 — 1073), it is difficult to be tested at LHC. We need pp colliders with 
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larger \/s. For example, if \/s = 100 TeV, with 3 ab™* luminosity, we can obtain O(10?—10*) 
events for pp > Z’ — hiha process in the mass region mz ~ (5 — 10) TeV which means 
it may become testable. CP-violation in the scalar sector will be confirmed if both nonzero 
hyZZ and Z'hıhə vertices are discovered. 

This model is attractive both theoretically and phenomenologically. Theoretically, in 
this model, we proposes a new possible CP-violation source, which may provide new under- 
standing of the matter-antimatter asymmetry problem in the Universe. Besides this, the 
spontaneous CP-violation mechanism is also a possible solution to the strong-CP problem, 
which is worthy to study further. This model is also a candidate to connect between the 
composite Higgs mechanism and CP-violation in the scalar sector. Based on this, new CP- 
violation effects are naturally suppressed by the global symmetry breaking scale f, as shown 
in the calculation of electron and neutron EDM. 

Phenomenologically, it is an application of the basic idea to measure hV V and Vhyho 
vertices. It provides an example to show how extra scalars and gauge bosons can help to 
confirm new CP-violation sources, which also implies the importance to search for VV S- 


and V SS-type vertices. It also shows another motivation for future ete~ and pp colliders. 
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Appendix A: Improved Formalism of the SLH Model 


In this section show the improved formalism for the SLH model based on [34]. The 
neutral scalar sector (including six degrees of freedom) can be divided into CP-even and 
CP-odd parts. The CP-odd part, denoting as G; running over 7, G, G’, and y?, is not 


canonically-normalized. We can write the kinetic term as 


1 
LD z Ku" Gð G3. (A.1) 
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The matrix elements of K are calculated to O ((v/f)*) in B4. If we rewrite this term in 


another basis S; = U;,;G; which is canonically-normalized 
1 p 
LJ 3 ou? SiOn Sj (A.2) 


thus we can define a inner product (5;|S;) = 6;; in the linear space spanned by the scalars 


S;. A straightforward calculation shows that 


(GiIG;) = (K); (A.3) 


The VEVs in ©; will lead to two-point transitions between gauge bosons and pseudo- 


scalars as 


LoVe F piOn Gi (A.4) 


P 


where V, denotes a gauge boson running over Z, Z’, and Y°, and F is a 4 x 3 matrix. The 


matrix elements of F are also calculated to O ((v/f)’) in B4. The gauge fixing term must 


provide the two-point transition like 


Lar. 2) (Ə V”) Fis (A.5) 


to cancel all contributions from (A.4). Define 


Gp = ere (A.6) 
in the convention of (which is also the convention of this paper), we can derive that 


(n|Gp) =0, and (Gp|Gq) = (Mọ) , (A.7) 


pq 


through a straightforward calculation where Mj, is the mass matrix for gauge bosons in the 


basis (Z, Z’, Y). Calculate to the leading order of (v/f) for every matrix element, we have 


v2 cow? v3 
Achy Ach, 3-12, 3V2ew tog f 
2 2f? v3 
M?, = g? caw Y —— l A.8 
V5I | Ig AB, 3—thy 3y/6—-12,,c2, toa f (A.8) 
63 u f2 
3V2cw top f 34/6—-13,,c2, toe f 2 


Using an orthogonal matrix 


R, we can diagonalize Mj, as 


RMy,R 
2 RT 


tian (7 —] 
pa Z pas and V, =|] 


Rpa Va, 


(A.9) 
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where V, denotes the mass eigenstate of a gauge boson and m, is its mass. The matrix 


elements of R are calculated to O ((v/f)*) in [84] as well. For simplify, to this order, the 


off-diagonal elements can also be expressed as 


(MB), 
Roe WE), = (MY, (A.10) 


It is natural for us to define 


A RygGq _ (R a Gi 


m 


(A.11) 


P 


According to (n|n) = (K7'),, © 1 + (2/t3,)(v/f)?, we should also define 


p 


n= W (A.12) 


It is easy to check that in the basis (i, Gr), the kinetic part is canonically-normalized. (A.5) 


also becomes mp (0,7) G thus it is natural to choose the gauge fixing term as 
1 m ~ \2 
Lar =~ zg (0 - mG) (A.13) 
p 


It is now clear that Ĝ, is the corresponding Goldstone eaten by V, and its mass should 
be V Elta where ép is the corresponding gauge parameter. To this step, we have already 
built the formalism to treat a model with non-canonically-normalized scalar sector and the 
SLH model is one of the examples. The main point is that all the two-point transitions 
must be carefully canceled if we don’t want these kind of Feynman diagrams appear during 
calculation. 


Because of the 7 components in the Goldstone fields, the interactions including 7 must 


be changed comparing with the naively calculated case. We divide F into 


F = @ *) (A.14) 


where fp = F, is a 1 x 3 vector and F is a 3 x 3 matrix. Thus for any kind of couplings 
including the pseudo-scalar degrees of freedom, if we write the coefficients as (cn, cj) in G; 
basis where c; runs for the couplings including G, G’ and y?, the physical coupling should 


be 


gea ie), (a =i (=) ) (A.15) 
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For example, the anti-symmetric type Vhn couplings in mass eigenstates can be parameter- 
ized as 
Levin = $ (dH — OD) (Cnu + Chin Zi, + an Ye) - (A.16) 
With the improved formalism, we can calculate to the leading order of (v/f) as 
gs 1 aa 272 (%\ eas 
“Eh 9/3 top (5) ad ee eee (5) ae ean) 


The first two results are quite different from those appearing in previous papers B5]. 


Similarly, the Yukawa couplings between 7 and SM fermions can be parameterized as 
imy -+ 
Lith =~ Ss ent dv fa. (A.18) 
f 


According to (A.15), cn, = 0 to all order of (v/f) for f = v,£,u,c,b. This result is also 
quite different from that in previous papers [35]. For f = t,d,s, we have 


=o v C28 + C269 Al 
Cn,f = —Of (=) a (A.19) 


which is generated by the left-handed mixing between SM fermion and additional heavy 


fermion. Formally all these results can be calculated to all order of (v/f), though some of 


the results are extremely lengthy. 
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